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Recently a new state of the bb system has been observed by the ATLAS Collaboration, with a 
mass of 10530±5±9 MeV. This state has been identified with the n — 3 P-wave radial excitations of 
£Sj . the bb system with parallel quark spins, called Xb(3P). The measured value of the mass corresponds 

to the average of the J = 0, 1 and 2 states, while the splittings of these states are not yet resolved. 
In this work we present predictions from different potentials models for the values of these splittings. 
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I. INTRODUCTION 



Heavy quarkonia, or mesons formed as heavy quark- antiquark bound states, constitute a valuable ground to study 
strong interactions, because they are related to both the perturbative and non-perturbative regimes of QCD. In heavy 
'■^ | quarkonia energies are marginally high enough for perturbative methods to become useful, while the confining, non 
perturbative aspect of the interaction, is also of crucial importance. 



O-i Recently a new state of the bb system has been observed by the ATLAS Collaboration, with a spin- weighted average 
( "barycenter" ) mass of 10530 ±5 ±9 MeV [lj. This state has been identified with the n = 3 radial excitation of the bb 
' 1 system with angular momentum numbers s = 1, £ = 1 (hence J = 0, 1, 2), called Xb(3Pj). The three states J = 0, 1, 2 
are closely spaced in mass, so the measured value of the mass corresponds to their average, while the splittings of 
| these states are not yet resolved. In this work we present predictions from different potentials models for the values 
■ of these splittings. 

Since the discovery of the J ftp meson in 1974, understood as a (cc) system, many potential models have been 
t-H ■ proposed to describe the interaction between the quark and antiquark in a bound state [2-25]. One of the first 
proposed models was the so called Cornell potential @, [H, @], which is a Coulomb-plus- linear term combination 
OO ■ that takes into account general properties expected from the interquark interaction: a Coulombic behavior at short 
distances and a linear confining term at long distances, representing the perturbative one-gluon exchange and the 
non-perturbative chromoelectric flux tube of confinement, respectively. Other models based on phenomenological 
grounds include a logarithmic potential [J and a non-integer power law potential [§]. Subsequently, elements from 
QCD were included in different ways in the potential formulation, from the inclusion of the running of the QCD 
coupling constant in the Coulombic interaction to a derivation of the short distance quark-antiquark potential 
from perturbative QCD, where spin-dependent interactions naturally appear [cl [Ici - fl5l ] . A modification of the Cornell 
potential to take into account a saturation effect in the linearly growing confining part has also been proposed, inspired 
in lattice data [TH, [ID, [2l|. In the context of these potential models, a successful description of many properties of 
quarkonia has been done, both for mass spectrum and decay widths. A comprehensive review of the status of heavy 
quarkonium can be found in Ref. |26| . 

Even though the b quarks are heavy compared to the scale of confinement in QCD, the b-b interaction is sufficiently 
strong to make the bound state somewhat relativistic. Quark model calculations give typical velocities v q ~ 0.25 — 
0.4. Consequently, one expects relativistic effects to be important, and at a fundamental level the very use of a 
potential, which represents an instantaneous static interaction, is questionable. Nevertheless, potential models have 
been successful in predicting mass spectra and decays in quarkonia, and are particularly useful in describing the mass 
splitting among different states of orbital angular momentum and spin. 

Here we study the prediction on the mass splitting in Xb{3P) from a set of potential models that have been able 
to predict other quarkonium spectra with relative success. We compare the average masses of known parts of the 
spectra in P-wave states with the results given by each of the potential models. An important issue to consider in the 
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model is the separation between scalar vs. vector nature of the quark-antiquark interaction, as the spin-dependent 
interaction is sensitive to this separation. Since this separation is purely phenomenological, we must first use as input 
a known part of the spectrum in order to fit the parameters that model this separation; we then obtain results for 
other parts of the known spectrum, which we use to check the validity of the model. Then, taking into consideration 
the degree of success of these results, we adventure in predicting the yet unknown mass splitting in the \b(3P) state 
within each of these models. 

In Section [TT] we briefly describe the potential models that we consider in our analysis. In Section Mil we present 
the results obtained from the different models including the predictions for the splittings of the Xb(3P) state. Finally 
in Section HVl we state our conclusions. 



II. POTENTIAL MODELS 

Rigorously speaking, a potential in a bound state is a static instantaneous interaction. Such a description is good 
only in the non relativistic limit. A systematic relativistic expansion of the interaction provides additional terms that 
are corrections that depend on the orbital angular momentum and spin of the bound particles. The expansion should 
represent charmonium and bottomonium rather well since the motion of their constituents is not highly relativistic. 
In what follows we consider four representative models given in terms of a spin- independent, static central potential, 
where relativistic corrections are incorporated as the aforementioned expansion. We call these models a) the Cornell 
model b) the Screened Cornell model c) the Richardson model, and d) the p-QCD (perturbative QCD) model, which 
we describe in what follows. In addition, for the first two models we consider two variations, henceforth called Cornell- 
I and Cornell-II, and Screened Cornell-I and Screened Cornell-II, respectively, depending on the separation of the 
potential into a vector exchange and a scalar exchange, as explained below. 

A. The Cornell model 

The Cornell potential 0] was originally proposed to describe masses and decay widths of charmonium states. It 
takes into account the following features of the interaction between the constituents: (i) a Coulomb-type interaction 
which should describe the short distance regime mediated by a single gluon exchange in perturbation theory, and (ii) 
a quark confining interaction that dominates at long distances. Explicitly, the potential is given by 

V(r) = -P + ^ + V , (1) 
r a 

where the first term corresponds to the Coulombic interaction and the second is the interaction responsible for 
confinement that arises from the color field flux tube between the quarks. The coefficients a and j3 are adjusted to 
fit the charmonium spectrum, but with the assumption that, roughly speaking, it should be valid for all other heavy 
quarkonia. As such, the flavor dependence should arise solely from the mass of the bound quarks. 

Many fits have been done to adjust these parameters [1,0, Hi}. In this work we will consider the most recent values: 
(3 = 0.52, a = 2.34 GeV -1 , m c = 1.84 GeV and mj, = 5.17 GeV 0, where /3, a and m c were fitted using as inputs the 
mass differences of the charmonium states J/^(15), ip{2S) and the Xc(l-P) spin-weighted barycenter, as well as the 
requirement that the charm quark mass should be close to the lightest charmed-meson mass. The value of mj, was 
chosen taking into account limits based on the T masses. 

The potential in Eq. (fT]) represents the static limit of the interaction. A relativistic expansion will include additional 
terms, that will constitute smaller corrections in the case of heavy quarkonia, and will depend on the orbital angular 
momentum and spin of the constituents. These terms are essential in order to reproduce or predict the mass splittings 
of quarkonium states that differ in angular momentum, such as the Xb(3Pj), for J = 0, 1, 2. 

A general expression for the spin-dependent potential of a qq bound state in QCD using a Wilson loop approach 
was calculated to order 1/m 2 by Eichten and Feinberg: [28[ 



V SD (r) = 



Si • Li S2 • L2 



2m\ 2m| 



dV(r) , 2 _dV^(r) 



rdr rdr 



(S 2 -Li-Si-L 2 )rfVS(r) 
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6mim 2 Am\mi 

where V(r) is the spin-independent part of the potential, while Vi(r), Vz(r) and V^ir) are the spin-dependent parts, 
that can be expressed as expectation values of color electric and magnetic fields. These potentials are not all inde- 
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pendent, as pointed out by Gromes [20|. Lorentz invariance imposes the relation V(r) + V\{r) 
eliminate Vi(r) from Eq. ([2]) and, in the equal- mass case, toi = m% = m, rewrite it as: 



V2(r). We can thus 



V SD (r) 
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2m 2 



dV(r) dV 2 (r) 
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(2S-S-3) V 2 y 2 (r), 



(3) 



where S = Sx + S 2 and L = = — L 2 . The terms in Eq. ([3]) are referred to as spin-orbit, tensor and hypcrfinc 
interaction, respectively. Notice that only the spin-orbit and tensor interactions are responsible for the mass splittings 
of the Xb(nPj). 

Now, the otherwise independent potentials V(r), V^r) and V^(r) are further related if one assumes that the quark- 
antiquark potential due to the strong interaction arises from the effective exchange of a scalar and a vector particle. 
This assumption is clearly part of the modeling, as there could be further effective exchanges, but if we stick to it, 
the quark-antiquark lagrangian has the form: 

Lint — s(q 2 )uuvv + v(q 2 )u^/ fJ ,uvj fJ 'v. (4) 

The expansion of this expression in powers of v 2 /c 2 gives a Breit-Fermi form of the spin-dependent potential [30l l3lj] : 

d 2 v(r) 
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(5) 



Comparing Eqs. ([3]) and ([5]), one can make the following identifications, which relate V(r), V2(r) and V^,(r): 

V(r) = v(r) + s(r), V 2 (r) = v(r), V 3 (r) = ^ - 

rdr dr z 

We then see that the spin-independent potential V(r) is composed by the vector and scalar parts, while the spin- 
dependent parts, V2<V) and V^{r), are given in terms of the vector part only. 

To leading order in the wave function obtained from the spin-independent potential V(r) (e.g. Eq. [TJ), one can 
express the spectrum of the Xb{ n Pj), for a fixed radial excitation n, as: 

M(P 2 ) = M + a- 26/5, 
M(Pi) = M-a + 2b, 
M(P ) = M-2a-46, 

where a and b are the expectation values of the radial functions in the spin-orbit and tensor terms, respectively: 

1 



(G) 



2m 2 
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12m 2 



ds ^ dv \ 
rdr rdr / ' 

dv d 2 v s 
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As such, the potential V(r) provides an explicit form for the sum v(r) + s(r), and for the bound state wave function to 
leading order. However, to obtain the explicit form of the spin-dependent potential, Eq. (j5j), it is necessary to know, 
or to define within the model, the separation of V(r) into the scalar and vector parts s(r) and v(r). 

For the Cornell potential, Eq. (TTJ, one can argue that a reasonable separation is to assign the Coulombic term to 
the vector part, coming from a single gluon exchange at short distances, and the linear term to the scalar part, coming 
from the flux tube at long distances: 



v[r) 



s(r) 



(7) 



Here we call this prescription the Cornell-I model. 

Alternatively one can relax this assumption by introducing two phenomenological parameters, r\ and A, that define 
the combinations of Coulomb and linear terms of the scalar and vector potentials (where the case r\ = A = 1 
corresponds to Eq. [7]): 



s(r) 
v(r) 



-(1-A)£+tj4, 

r 
r a z 



(8) 
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One can then find the values of rj and A that provide the best fit to a part of the spectrum. In this study we look for 
the values that minimize the following % 2 function, using the mass splittings in the experimentally known XbO-Pj) 
andxfc(2Pj): 

2 _ v ( AM ij -(rcP)( modcl > - AM^jnP)™ \ 2 
X k, \ ^(nP)(ex P ) j > W 

y =21,10 

where 

AMai(lP) = M(LP 2 ) -M(lPi), AM M (1P) s M(li\) - M(1P ) (10) 

are the mass splittings in the Xi,(lPj), and similarly are the AMy(2P) defined for the Xb{2Pj)- Following this 
procedure, we find the optimal values rj = 1 and A = 0.923. We call this prescription the Cornell-II model. 

It is interesting to notice that the optimal values of rj and A are remarkably close to the Cornell-I prescription, 
which is the one commonly assumed in the literature. 



B. The Screened Cornell model 



A variation of the Cornell potential, which we call the Screened Cornell potential, has been used to include the 
effect of saturation of the strong interaction at long distances that appears in lattice data jig : 

\ r cry \ fir J 

where fi is the screening parameter. This potential behaves like a Coulomb potential at short distances but, unlike 
in the previous model, it tends to a constant value for large r (namely, for r 3> fi 1 )- In other words, the linearly 
growing confining potential flattens to a finite value at large distances, corresponding to the saturation of as to a 
finite value for decreasing Q 2 (321-34] . This effect should be due to the creation of virtual light quark pairs that screen 
the interaction between the bound quarks at long distances. The values of j3, a and fi we use here are those of Ref. 



21 [ and are shown in Table HI Of these values, (3 and /j, are intrinsic to the model, while a, m c and nib were fixed by 
the authors in order to reproduce the J/ip mass in cc and the T(1S) and T(2S) in bb. 

For the potential in Eq. (jlip , we can separate the scalar and vector parts in the same common way as in the Cornell 
potential (that is, the Coulombic term as vector and the linear term as scalar), a prescription we call here the Screened 
Cornell-I model. 

Altenatively, just as before, we can introduce two phenomenological separation parameters, r\ and A, which are 
fitted by the same optimization procedure: 

*)- (-d-^ + ,4)(^ 

V r a z J\fir 

v(r) = + (1 - r;)^) (^—) ■ (12) 

\ r a z / \ \ir J 

The optimal values we find in this prescription are 77 = 0.810 and A = 1. We call this prescription the Screened 
Cornell-II model. 



C. The Richardson model 



The Richardson potential is another well known model that incorporates the features of asymptotic freedom at short 
distances and linear confinement at long distances @. With a minimal interpolation between these two asymptotic 
behaviors, the Richardson potential is obtained: 



V(r) 



8tt 
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-A 



Ar 



/(Ar) ' 
Ar 



(13) 
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Table I: Best values of the parameters for each of the potential models, according to previous authors, except mb for the 
Richardson model, which was updated with the current experimental data. 





Equation 


P 


a GeV -1 


m b GeV 


V GeV 


fi fm 1 


A GeV 


Cornell 


© 


0.52 


2.34 


5.17 


-0.50805 






Screened Cornell 




0.423 


1.858 


4.6645 




0.71 




Richardson 


G3 






4.8877 






0.398 



with 

/(t)=1 . 4 r^ . 

W A 3 ln 2 (g 2 -l)+^ 2 

Here n/ is the number of light quarks relevant to the renormalization scale, taken equal to three, while A = 0.398 GeV 
is the scale of interpolation between the two asymptotic regimes and m c = 1.49 GeV is the constituent mass of the 
charm quark in the model. These values have been fitted to reproduce the mass of the charmonium states J ftp (IS) 
and ip(2S). For bottomonium the same value of n/ and A should be used, setting the bottom quark mass at 
rrih = 4.8877 GeV in order to reproduce the current T(15) mass. According to Richardson, the value of nj is kept 
equal to three (i.e. the number of light quarks), since the Appelquist Carazzone theorem [35| implies that the effect 
of quarks heavier than the energy scale that determines the dynamics should be small (the latter is related to the 
binding energy or the inverse of the radius). With these values for A, rif and mb, the T(2S) mass is predicted with 
good agreement with experiment. 

In contrast to the Cornell potential, the Richardson potential has no obvious separation into scalar and vector 
parts. In order to be able to determine the splittings with this potential, one needs this separation. The following 
phenomenological separation has been used in previous works |31| : 

v(r) = V(r)e' r2/a ', 

s(r) = V(r)(l-e- r2/a '). (14) 

Here V(r) is the Richardson potential of Eq. (fT3"]) and a r is a phenomenological length scale that separates the vector 
character at short distances and scalar character at large distances. Following the same % 2 optimization procedure as 
before, the value a r = 0.1782 fm is found. 



D. The p-QCD model 



One last potential model we consider in our study (121 136[ is based on a semi-relativistic treatment of perturbative 
QCD interactions to one loop, which we call the p-QCD model, where the spin-dependent potentials, which appear 
in the perturbative treatment, are somewhat different than those in the previous models. The Hamiltonian in this 
model considers a relativistic expression for the quark kinetic energy, Hq — \J p 2 + m 2 , unlike the previous models 
where a purely non-relativistic approximation, Hq = m + p 2 /2m, is used. Additionally, the potential is composed of a 
short distance part obtained from perturbative QCD to one loop, and a phenomenological long range confining part. 
The complete potential has the form 



1- ^ + ^(33-2n/) (In M r + 7E ) 



Ar + V S + V L . 



(15) 



The first term is the spin independent Coulomb-like part, corrected to one loop in QCD. The term Ar, where A is 
a phenomenological constant, is the spin-independent, long range confining potential, assumed to be linear in r. Vs 
is a short-range part of the potential that appears in the perturbative calculation to one loop, and includes most 
spin-dependent parts in the form of spin-orbit (LS), tensor (T) and hyperfme (HF) terms, and a short-interaction 
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spin independent term (SI): 
V LS 



2a s L-S j as 
6tt 
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Finally Vt is the relativistic correction arising from the confining potential (thus proportional to A), also composed 
of spin-dependent terms: 



V L = -(l-fv) 



A 



2m 2 r 



L • S + h 



2m 2 r\^t^ 



JL4_ 

2m 2 i 



L S 



A 



3m 2 r 



(3Si f S 2 • f - Si • S a ) 



The phenomenological parameter fy was introduced by the authors to represent the fraction of vector vs. scalar 
character of the confining potential, which they fit with the available data. For bottomonium t he p arameters A, as, 
to, fi and fy were fitted using as input eight masses of the bb spectrum. The authors obtained [33] A = 0.175 GeV 2 , 
a s = 0.295, to = 5.33 GeV, fi = 4.82 GeV and f v = 0. 

Remarkably enough, fy = implies, just like in the other models, that the long distance part of the potential 
corresponds to purely scalar exchange. 



III. CALCULATIONS AND RESULTS 



A. The spin-average masses of \b(!P), x>(2P) and Xb(3P) 



As seen in Eq. ([6]), the masses of the states Xb(nPo), Xb(nPi) and XbinP^) differ by small amounts. These differences 
are determined by spin-dependent interactions. On the other hand, the barycenter (or spin-averaged) mass of these 
states is defined as 



M 



5M(P 2 ) + 3M(P 1 ) + M{P ) 
9 



(17) 



and is determined in each model by the leading, spin-independent part of the potential. Consequently, we expect 
the spin- independent part of the models we previously described to correctly reproduce the experimental baryecntcrs. 
We have calculated M for the Xb0-P)> Xb(2P) and Xf,(3P), numerically solving the Schrodinger equation 38( for the 
potential models described in the previous section, using the parameters listed in Table U for the first three models, 
respectively. The results we obtain are shown in Table [TTJ We include the results obtained in Ref. [37j by the authors 
of the p-QCD model. The value of the constant Vq for the Cornell potential [Eq. ((TJ] was fitted in order to obtain the 
experimental value of the T(IS'). 

We should recall that the parameters in the Cornell and Richardson potentials, which are shown in Table U were 
fitted using charmonium states (except for the bottom constituent mass rrib), and those of the Screened Cornell model 
were fixed using both cc and bb masses. In contrast, in the p-QCD model the free parameters were fitted using bb states 
only. Consequently, one should naturally expect the latter model to give a closer prediction for the Xb, in comparison 
with the other models, where an accurate prediction for Xb masses is clearly a more demanding requirement. 

The bottom line of Table [TT] shows the experimental values of the barycenter masses, in order to compare them with 
the model results. To ease the comparison we also include, for each Xb, the mass discrepancy 5 n p between the model 
result and the corresponding experimental value. 

As a first observation from Table [TT1 one can see that in most cases the model discrepancies are larger than the 
experimental uncertainties, indicating that in general there is still need for improvement in the models. 
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Table II: The mass barycenters of the states Xb{lP), X&(2P) an d Xi>(3P) reproduced by the potential models described in 
Section |TIJ and compared to the corresponding experimental values. On the right of each reproduced mass is the discrepancy 
between this value and the corresponding experimental value, 5 n p = M(nP) ( moc ) e i) — M(nP)( ex py The experimental values 
for IP and IP are from Ref. [39l ]. while those for 3P are from Ref. [l[. All values are in MeV. 



Model M(1P) 


Sip 


M(2P) 


&2P 


M(3P) 


83P 


Cornell 9958.3 


58.4 


10312.6 


52.5 


10595.3 


65.3 


Screened Cornell 9907.9 


8.0 


10261.2 


1.0 


10516.4 


-13.6 


Richardson 9895.7 


-4.2 


10248.8 


-11.4 


10520.1 


-9.9 


p-QCD [37J 9898.7 


-1.2 


10261.2 


1.0 


10543.9 


13.9 


Experiment 9899.87 ± 0.27 




10260.20 ± 0.36 




10530 ± 10 





Notwithstanding, the errors are quite small in all models except Cornell. This tendency indicates that the inclusion 
of saturation (i.e. going from Cornell to Screened Cornell) clearly improves the predictions for the barycenter masses, 
but there is no clear distinction in the reliability between the last three models at this stage. 

We must also comment that the discrepancy in the Cornell model can be considerably improved in a very simple 
way: the offset Vq was fitted using the T(IS') mass. If we had fitted Vb using \b(lP) instead, the discrepancy in the 
barycenter mass of Xb(?P) and Xb(3P) in Table [TT] would have been 5 2 p = -6.1 MeV and S 3P = 6.8 MeV, instead 
of 52.5 and 65.3 MeV, respectively. In this sense, the Cornell potential does not fit very well the whole bb spectrum, 
but could still be good in predicting mass differences within Xb{nP) states, which are all of the same spin and orbital 
angular momentum. 



B. The splittings of Xb(lPj), Xb(2Pj) and Xb(SPj) 



To date, the splittings of the J = 0, 1, 2 states in the Xb(lPj) and Xb(2Pj) have been experimentally measured with 
reasonable precision, while there is still no data for the splitting in the Xb(3Pj)- Our main goal here is to predict the 
latter splitting according to each of the aforementioned models. 

We characterize the splittings of the J = 0,1,2 states within a given Xb(nP) by the mass differences AM2i(nP) 
and AMio(nP), as defined in Eq. (flU)) . An additional parameter which is used in the literature to characterize the 
splittings is the ratio of the mass differences within a given Xb{nP), defined as 



AM 21 (nP) 
AMi (nP)' 



R x (nP) = (18) 



In our calculations, we first reproduce within each model the splittings in the Xb{^Pj) and xt(2Pj) and compare 
these results with their corresponding experimental values, as a way to check the reliability of the models. These 
results are shown in Tables UlTl and HVl respectively. In these tables one can make several observations. 

First, from Table HTT1 for the splitting in Xb(lP), the best results for the AMy are obtained with the Richardson and 
p-QCD models, followed by the Cornell-II model, while the largest discrepancies appear with the Screened Cornell-I 
model. On the other hand, concerning the splitting ratios R x , the pattern of performance is different: the Screened 
Cornell-I joins the Richardson and p-QCD models in giving the best results, while the largest discrepancy occurs with 
the Screened Cornell-II model. 

Alternatively, we can see whether the results for the splitting in Xb(lP) improve in going from the Cornell to 
the Screened Cornell model (i.e. by including saturation at long distances), as it happened in the results for the 
barycenters. Clearly this is not the case: there is no such an improvement on the splittings: the Screened Cornell 
models actually perform worse than the simple Cornell models. 

Finally, we can see whether the best fitted vector vs. scalar separation causes an improvement in the results (i.e. 
going from models of type -I to type -II) . Indeed, a moderate improvement can be seen in Table IIIII going from 
Cornell-I to -II or from Screened Cornell-I to -II, which in any case is expected, because the fitting was done with 
that purpose. However, the improvement occurs only in the AMy, but not so in the ratio R x . 

In turn, for the splittings in x;,(2P), the best results for AM y - are given by the p-QCD model, and Screened Cornell- 
I and -II models, while the largest discrepancy is found in the results given by the Cornell-I model. In contrast, for 
the splitting ratio R x the best results are given by the Screened Cornell I and p-QCD models, while the largest 
discrepancy occurs in the Screened Cornell-II model. 
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Table III: The mass splittings in the Xi,(lP) mesons calculated for each of the models described in Section [TTJ and their 
experimental values: the mass differences AM21 and AM10 are defined as in Eq. (|fOp. and the splitting ratio R x is defined as 
in Eq. (pi). 

Model AM 2 i(IP) [McV] AM 10 (1P) [MeV] R X (IP) 



Cornell-1 23.13 32.24 0.717 

Cornell-II 20.12 29.14 0.690 

Screened Cornell-1 16.19 25. 2f 0.642 

Screened Cornell-II 20.87 28.46 0.733 

Richardson 18.58 30.87 0.602 

p-QCD [37] 19.3 29.94 0.645 



Experiment 19.43 ± 0.57 33.34 ± 0.66 0.583 ± 0.021 



From the above one sees that the results are rather disperse: models that reproduce well some of the splitting 
features do not do so well on other features. The inclusion of saturation in the Cornell model does not seem to 
make it much better, and the optimal adjustment of the vector vs. scalar separation does not cause a significant 
improvement either. 

The sole exception, to some degree, is the p-QCD model, which tends to be more often within the best results. 
However, this tendency was expected, since this model contains a larger number of parameters and they were fitted 
purely with bb. In contrast, all the other models were built using the charmonium spectrum and are now required to 
fit the bb spectrum as well. 

One clear pattern in the results for Xf,(lP) and x&(2P) is that all models give values for the ratio R x which are 
consistently larger than experiment. This feature could be an indication that the separation of the potentials purely 
in terms of scalar and vector exchanges may not be sufficient, and an expansion into exchanges of additional spin and 
parity may be required [401 ] - 

Now, concerning the prediction of the models for the so far unknown splitting in the Xb(3-P./), Table fVl shows the 
results of our calculations, together with the prediction of the p-QCD model according to its authors [4l|. 

Looking at the values for the ratio R X (3P), if we extrapolate from the results on Xb(l-P) and Xb(2P) we may suspect 
that these are also overestimations, i.e. all models tend to give too large values for R x . Within that assumption, we 
may consider more successful those models that give the smallest prediction for R Xl which correspond to the p-QCD 
and the Screened Cornell-I model. 

However, if we look at the AA/y results in the 3P, the predictions of the Screened Cornell-I model are much lower 
than all the others. Then, going back to IP and 2P, where there is experimental evidence, this model also gives too 
small splittings in those states. In that sense, concerning the splittings the reliability of this model is questionable. 

Besides that observation, among the other models there is a tendency in AM10 to be near 22 MeV. For AM21 the 
predictions are more dispersed, but centered around a value of 15 MeV. If we consider these two values as the best 
predictors for the splitting, one then predicts a ratio R x = 0.68 for the x&(3P). 

Using the definition of the barycenter given in Eq. (TTTl) we can deduce the masses of the states J = 0, 1 and 2 of 
the X(3Pj) in terms of the barycenter and splittings: 



Af(P ) = M — 



5AM 21 + 8AM ; 



10 



9 



M(P = m- 5AM " 9 - AM "> , (19) 

M(P 2 ) = M+ 4AM21 Q +AMl ° . (20) 
9 

Using the experimental value and its uncertainty for the barycenter and the estimated splittings, the masses of the 
three x 6 (3P) states (J = 0, 1, 2) would then be 10502 ± 10 MeV, 10524 ± 10 MeV and 10539 ± 10 MeV, respectively. 
On the other hand, the threshold for B + B~ decay is 10558.5 ± 0.3 MeV. Consequently our estimate of the J = 2 
state is just 20 ± 10 MeV below the threshold. 
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Table IV: The mass splittings in the Xb{2P) mesons calculated for each of the models described in Section [TTJ and their 
experimental values: the mass differences AM21 and AM10 are defined as in Eq. (JTDJ, and the splitting ratio R x is defined as 
in Eq. tfTSJ. 



Model AM 2 i(2P) [MeV] AMwQP) [MeV] R X (2P) 

Cornell-I 19.22 26.43 0.727 

Cornell-II 16.77 23.91 0.701 

Screened Cornell-I 12.91 19.34 0.668 

Screened Cornell-II 15.91 21.45 0.742 

Richardson 18.27 25.36 0.720 

p-QCD [37] 16.40 24.30 0.675 

Experiment 13.5 ± 0.6 23.5 ± 1.0 0.574 ± 0.035 



Table V: The predictions for the mass splittings in the Xb(3P) mesons according to each of the models described in Section HT1 
the mass differences AM21 and AM10 are defined as in Eq. (|10[) . and the splitting ratio R x is defined as in Eq. (|18|) , 



Model AA/ 2 i(3P) [MeV] AMi (3P) [MeV] R X (3P) 

Cornell-I 17.7 24.0 0.734 

Cornell-II 15.4 21.8 0.710 

Screened Cornell-I 10.8 15.8 0.683 

Screened Cornell-II 12.9 17.3 0.747 

Richardson 17.1 22.1 0.773 

p-QCD [41] 148 22T 0.670 



IV. SUMMARY AND CONCLUSIONS 



Recently, the Xb(3P) state was observed and its barycenter was determined, while there is still no available data 
regarding its mass splittings. 

Here we have used four well known potential models to predict the so far unresolved splitting in the Xb(3P) states. 
The models were described in Section ITTI and are called here the Cornell, Screened Cornell, Richardson and p-QCD 
models, respectively. The first three models were originally adjusted to reproduce the charmonium spectrum, and are 
required to work for bottomonium as well. In contrast, the latter model was fitted using bottomonium states only. 

As a first attempt to test the reliability of these models, we use them to reproduce the known barycenter masses of 
the Xb{lP), Xb(2P) and Xb(3P) states. 

The Cornell model, which is the simplest of all of them, shows the largest deviations in its results for the barycenter 
masses. The Screened Cornell model, which is similar to the previous one, but with the additional feature of saturation 
at long distances, shows a clear improvement in its results for the barycenters. The other two models, Richardson and 
p-QCD, give similar values. Up to here, it seems that the latter three models perform similar in terms of reliability. 
In any case, one should notice that the discrepancies in these model results are in general larger than their current 
experimental uncertainties. This is an indication that there is still need for model improvement. 

Now, concerning the splittings, we recall that the barycenters of the Xb states are determined by the spin-independent 
part of the potential, while the splittings are determined by the spin-dependent part. The latter, which is a relativistic 
correction, should be sub-dominant for heavy quarkonia like the Xb- The spin-dependent part, however, is not 
unambiguously determined as an expansion of the potential, as it also depends on the Lorentz character (vector vs. 
scalar) of the effective interaction. In the models we study here, we have used a common phenomenological treatment 
of this feature, which is a separation of the potential into a vector and a scalar exchange, where it is usually assumed 
that the vector part should dominate the short distance regime, while the scalar part should represent the confining 
regime at longer distances. 

Accordingly, in the cases of the Cornell and the Screened Cornell models, we considered two variants for the scalar 
vs. vector separation. Our variant I considers the commonly assumed separation where the vector part corresponds 
to the short distance Coulombic term and the scalar part corresponds to the long distance linear (or linear-saturated) 
term. Our variant II uses an optimal scalar-vector separation that minimizes the error in reproducing the known 
splittings in the XbO-Pj) an d Xb{2Pj)- Interesting enough, we found that this best fitted separation is very close 
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to the usually assumed separation in variant I, in which the short distance part is purely vector exchange and the 
confining long distance part is purely scalar. This result also coincides with a similar analysis done in the p-QCD 
model by its authors, where they found that their best fit corresponds to a long distance part being purely scalar. 

For the Richardson potential, unlike the previous models, there is no obvious separation into scalar and vector 
parts, so we introduced a purely phenomenological function with one parameter to carry out a short distance vs. long 
distance separation. This parameter, also fitted with the Xfc(l-P) and Xb(2P) states, is a length scale that separates 
the two regions; the short distance part is then assumed to be mediated by vector exchange and the long distance 
part by scalar exchange. 

Of all the models studied, those that give more consistent results for the splittings are the p-QCD model and, to 
some degree, the Richardson model. One could have expected the p-QCD model to give the best results because 
of its more detailed short-distance treatment from perturbative QCD, and because it contains a larger number of 
parameters, all fitted to the bottomonium spectrum. 

Concerning the comparison between the Cornell and Screened Cornell models, we would have expected an im- 
provement in the reproduction of the splittings in the latter, as was the case for the reproduction of the barycenters. 
However it was not the case: we found the Screened Cornell model to give results that are no better than the simpler 
Cornell model. On the other hand, when comparing the variants I and II of these models, we had expected better 
results in the variants II, because they were fitted to do so. This was indeed the case in the results for the AMy, but 
not so for the ratios R x . 

The one consistent pattern in all results for the splittings is that all the models give ratios R x in Xb(lP) and Xf,(2P) 
that are larger than the corresponding experimental values. Concerning the reliability of the models, this tendency 
may be an indication that the composition of the potential in a scalar and a vector part may not be enough, and 
exchanges of further spins and parities may be present. From the point of view of predictions, on the other hand, 
we may then expect that the predictions for the Xb(3P) splittings will follow that tendency, i.e. the model with best 
predictions should give the smallest value for R x . This tendency again falls on the p-QCD model. 

Now, looking at the predictions for the AMio and AM21 in Xh(3P), one can see that the largest deviations around 
a central tendency are found in the Cornell-I and Screened Cornell-I models. Besides those predictions the tendency 
in AM10 is to be near 22 MeV, and for AM 2 i the predictions, although more dispersed, are centered around a value 
of 15 MeV. With these estimates and the experimental value for the barycenter, the masses of the three %f,(3P) states 
(J = 0, 1, 2) would be 10502 ± 10 MeV, 10524 ± 10 MeV and 10539 ± 10 MeV, respectively (the model uncertainty is 
not included). This means that the J = 2 state would be just 20 ± 10 MeV below the B + B~ threshold. 

Finally, if we consider these values for AM21 and AM10 as the best predictions for the splitting, one then predicts 
a ratio R x = 0.68 for the x&(3P), which is again quite close to the result of p-QCD. 
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